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ABSTRACT 

The establishment of a quasi-steady flow in a generic scramjet combustor is studied for 
the case of a time varying inflow to the combustor. Such transient flow is characteristic of 
the reflected- shock tunnel and expansion-tube test facilities. Several numerical simulations 
of hypervelocity flow through a straight-duct combustor with either a side-wall-step fuel 
injector or a centrally-located strut injector are presented. Comparisons are made between 
impulsively started but otherwise constant flow conditions (typical of the expansion-tube 
or tailored operation of the reflected-shock tunnel) and the relaxing flow produced by the 
“undertailored” operation of the reflected-shock tunnel. Generally the inviscid flow features, 
such as the shock pattern and pressure distribution, were unaffected by the time varying inlet 
conditions and approached steady state in approximately the times indicated by experimental 
correlations. However, viscous features, such as heat transfer and skin friction, were altered 
by the relaxing inlet flow conditions. 
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1. INTRODUCTION 


The National Aero-Space Plane (NASP) project has been the focus of the recent revival 
of hypersonic aerodynamics and propulsion research. Propulsion studies for the high Mach 
number range ( Mfn g h t > 4) have concentrated on the air-breathing supersonic-combustion 
ramjet (scramjet) using hydrogen as the primary fuel. Most of the early research studies 
of the mixing and combustion processes in a scramjet relied on conventional ground test 
facilities in which the test gas is heated to flight enthalpy by combustion of hydrogen or with 
an electric arc. These tests have generally been limited to Mfu g ht < 8 because of the limited 
total temperatures available with these facilities. Generation of test flows with enthalpies 
relevant to the flight regime, 10 < Mfu ght < 25, requires a different approach. 

One type of facility capable of this high enthalpy range is the pulse- type facility in which 
the test gas is rapidly heated by the passage of a strong shock wave. Two examples of 
pulse facilities are the reflected- shock tunnel and the expansion tube. A review of these 
facilities and their use in experimental hypervelocity aerodynamics has been presented by 
Stalker [1]. When using pulse facilities to simulate high enthalpy flows, a number of trade- 
offs between test time and test gas conditions must be made. One of the difficulties in using 
pulse facilities to test supersonic combustors is that the test times are relatively short - on 
the order of one-half millisecond. These short test times are a concern because, even in 
steady flow, many of the important flow processes in the combustor will take significant time 
to become established (i.e., approach steady state). With a transient flow the establishment 
may be further delayed. 

Much effort has been directed toward the optimization of the scramjet combustor and, 
in particular, to the efficiency of fuel-air mixing and reaction so as to obtain minimum 
length and wall heat load in the engine. Because of the high velocity, high temperature and 
low density of the combustor flow, the mixing of the fuel and air is relatively slow at the 
higher flight speeds. Consequently, mixing augmentation and enhancement by the controlled 
generation of turbulence and/or large-scale vortices in the mixing and flame-holding zones are 
being examined in both numerical studies [2] and experimental studies [3]. These techniques 
generally involve some type of recirculating flow in the wake of a body, step or transverse 
jet and hence require some finite time to reach steady state (or a quasi-periodic state). 

The focus of this study is to examine several time accurate simulations of steady and 
transient flow through a generic scramjet in an attempt to assess the following issues: 

• The time required for the flow in the wake region of the fuel injector (without injection) 
+ o become established relative to the typical test times available in pulse facilities. 
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• The impact of test flow relaxation (as encountered in the undertailored operation of a 
reflected-shock tunnel) on the approach to steady state combustor flow. 

• The validity of using the “hypersonic equivalence principle” (see, e.g., [4]) to transform 
the relaxing flow to a quasi-steady flow. 

The flow through the supersonic combustor was computed at the high enthalpy, short 
duration, and sometimes transient conditions typical of a shock tunnel or expansion tube. 
Calculations of the combustor flow were made in a time-accurate manner using the Navier- 
Stokes code SPARK [5], [6]. This code has been previously applied to the simulation of 
scramjet flows with fuel injection and a thrust ramp in a shock tunnel environment [7]. In 
the present study the code is used to compute the basic combustor flow without fuel injection. 

2. OPERATION OF PULSE FACILITIES 

A brief description of the operation of two particular pulse facilities will be given in this 
section. The first is a ffee-piston driven reflected-shock tunnel, T4, located at the University 
of Queensland [8]. The second is an expansion tube, HYPULSE, located at General Applied 
Science Laboratories, New York [9]. 

2.1. Reflected-Shock Tunnel 

The principal features of a shock tunnel and its operation are shown in Figure 1. The 
driver tube, which initially contains low pressure helium, and the shock tube which contains 
the test gas are separated by the primary diaphragm. Attached to the downstream end of 
the shocktube is a nozzle whose throat is significantly smaller than the diameter of the shock 
tube. The test gas is contained by a thin mylar diaphragm which separates the shock tube 
from the evacuated test region. 

The first stage of operation of the T4 facility is the launch of the piston and its acceleration 
down the compression tube. The driving force is supplied by the compressed air from the 
reservoir. The primary diaphragm, which is typically composed of two sheets of 2mm mild 
steel, subsequently bursts at a pressure 56.6 MPa. At this point, the helium has been 
compressed to 1 / 60th of its initial volume and is contained in approximately 0.5m of the 
compression tube. 

After diaphragm rupture, helium driver gas expands into the shock-tube and compresses 
the low pressure test gas before it. The incident shock travels the length of the shock 
tube, reflects from the closed end of the tube and brings the test gas to rest in the nozzle 
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supply region. Operation in this manner is called “tailored” [10] and is shown in the (x- 
t) diagram (Figure 1(b)) by the contact surface coming to rest when intercepted by the 
reflected shock. The compressed test gas is contained in the nozzle supply region in a length 
of approximately 0.25 — 0.5m. Ideally, the nozzle supply conditions are maintained as the 
reflected shock continues upstream through the driver gas. 

Upon shock reflection, the light secondary diaphragm bursts and some of the compressed 
test gas expands into the nozzle. Smith [1 1] has provided a quasi-one- dimensional model for 
the nozzle starting processes. 

For combustion studies, the useful test time is terminated by the contamination of the test 
gas by the driver gas. The mechanism for this contamination is the bifurcation of the reflected 
shock into two weaker oblique shocks near the wall of the tube. This mechanism has been 
studied in [12] and has been shown to occur for undertailored conditions at high enthalpies 
[13]. For undertailored operation, conditions are such that, when the reflected shock reaches 
the helium-air interface, it accelerates into the helium and an expansion propagates into the 
nozzle supply gas. With no other influences, the supply region will reach a new equilibrium 
but there will be a significant and unavoidable drop in nozzle supply pressure P t shortly 
after shock reflection. The expansion of the test gas slug delays the arrival of helium jetting 
down the walls of the shock tube and prolongs the useful test time. The approximate time 
at which contamination of the test gas is expected to reach the nozzle throat decreases with 
increasing stagnation enthalpy, H ly (see, e.g., [1]). 

The effect of finite driver size upon the relaxation of P t is also noticeable. Stalker [14] 
has suggested a mode of operation where the motion of the piston (after rupture) maintains 
approximately constant conditions in the driver tube while nearly 50% of the helium flows 
into the shock tube. Following this initial period, of approximately 1ms, the pressure in the 
driver decays rapidly and this decay is transmitted downstream along a (u+o) characteristic. 
The effect of the finite volume of the driver tube is that the the decay in P, may start early 
in the test period. 

The net result is that the history of the nozzle supply pressure and a representative pitot 
pressure at the exit plane of the Mach 5 nozzle will appear as shown in Figure 2. The 
pressure transducer measuring P, is located a few centimetres upstream of the reflecting end 
of the shock tube. Its trace shows the passage of the incident shock and the arrival of the 
reflected shock. Due to transducer rise time and location, the full reflected pressure is not 
recorded. Past the peak value, there is a continuous decay in P, due to the combined effects 
of undertailored operation and driver pressure decay. For helium driver gas and operation 
considered here, this decay is typically 25 — 30% during a nominal 0.5ms test time [15]. Note 
the delay between the supply pressure trace and the pitot pressure trace. This delay will be 
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used in the normalization procedure discussed in Section 3.5. 

2.2. Expansion Tube 

An expansion tube facility consisting of a driver tube (initially containing high pressure 
driver gas), a shock tube (containing the test gas) and an acceleration tube (containing low 
pressure acceleration gas) is shown schematically in Figure 3. The driver tube and shock 
tube are separated by the primary diaphragm and the acceleration tube is separated from 
the shock tube by a secondary (and very light) diaphragm. The shock- and acceleration 
tubes have the same diameter. Although the HYPULSE facility has a fixed driver tube 
rather than a free-piston driver, the gas-dynamic processes in the shock tube are similar. 

The operation is initiated by the rupture of the primary diaphragm. High pressure driver 
gas flows into the shock tube and the incident shock compresses and accelerates the test gas. 
On reaching the end of the shock tube, the secondary diaphragm bursts (ideally without 
causing disturbance) and the shock-compressed test gas expands into the the low pressure 
gas in the acceleration tube. The pressure to which the acceleration gas rises is below the 
pressure of the bulk of the shock-compressed test gas and so the downstream portion of the 
test gas undergoes an unsteady expansion to the test flow conditions. 

The test time commences after the test-gas/ acceleration-gas interface arrives and usually 
finishes with the arrival of the downstream end of the unsteady expansion fan. Figure 4 
shows a typical history of the static pressure in the test region. The trace shows a rapid rise 
due to the shock through the acceleration gas, a slower rise (possibly due to diffuser starting 
processes) to the test conditions at approximately 60/is and then roughly constant until the 
arrival of the unsteady expansion 400/ia later. 

3. NUMERICAL SIMULATION 

The computations of flow through the generic supersonic combustor were made with the 
Langley Research Center SPARK-2D code. This two-dimensional, elliptic, finite difference 
code was developed by J. P. Drummond [5], [6] to integrate the conservation equations for 
mass, momentum and energy in a time accurate manner. Since the primary concern is the 
temporal development of the combustor flow features, the scope of the modelling was reduced 
by considering only laminar flow of a nonreacting test gas and avoiding the fuel injection and 
mixing issues. The effects of transition, turbulence, mixing and chemistry will be examined 
in future studies. 
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3,1. Combustor Geometry and Boundary Conditions 


A generic scramjet combustor, typical of those tested in the reflected shock tunnel T4 
[15], is shown in Figure 5. This model consists of the major features of a scramjet combus- 
tor, including some form of injector (with an associated wake/ mixing region) and a set of 
confining walls (with their associated boundary layers). 

To perform the two-dimensional computations, the scramjet combustor duct in Figure 5 
was modelled as shown in Figure 6. Only half of the duct was considered in the computational 
domain. Two classes of calculations were performed which differed in the location of the 
inflow boundary. Cases 1,2,3 and 6 with the inflow boundary at the start of the duct walls 
(as = 0), have a flow domain as shown at the top of Figure 6. These are called the “short” 
injector cases. Cases 4 and 5, with the inflow boundary at the injector strut leading edge 
(x = -0.138m), are called the “long” injector cases (Figure 6(b)). Descriptions of the six 
cases considered in this study are shown in Table 1. 

Boundary conditions on the wall (at y = 0.02357 m) and the surfaces of the injector were 
set to zero velocity, zero normal pressure gradient and a fixed temperature T = 300JT (i.e., 
cold wall). Cases 4 and 5 included a free boundary along y = 0.02357m for —0.138m < x < 
0 m. Conditions at the inflow plane were supersonic. The outflow boundary at x = 0.5m 
was specified as supersonic with zero gradients in the flow direction. In Cases 1 and 2, the 
x-axis downstream of the injector (x > 0.114) was set as a cold wall so as to provide a 
situation close to that found in parallel wall injection studies. Cases 3,4,5 and 6 set this line 
to be a symmetry plane with zero normal velocity and zero normal gradients for pressure, 
temperature, and tangential velocity. 

3.2. Grids 

The grids used in the computation are shown in Figure 7. The entire domains, including 
the injectors, were included in the grids and during the computation all grid points within 
the injector were reset to initial conditions. Each grid consisted of a number of separate 
zones in which the nodes were clustered at one or more of the zonal edges so that large flow 
gradients near the boundaries and planes defining the injector could be resolved without 
excessive computational expense. 

The x and y coordinates for the grids used in the short injector cases were independently 
generated using one of Roberts’ [16] stretching transformations. (See also [17], Section 5-b.l). 
In each zone there were N + 1 nodes, including the end points, located at 

z = Z a fj + z b {\ - fj), ( 1 ) 


5 



where 


( 2 ) 


[(/? + 2a)A - /3 + 2 a] 
( 2 a + 1) [1 + A] 

(r,-a)/(l-a) 




V = jj,i = 0-N. 


( 3 ) 

( 4 ) 


Details of the zonal boundaries and stretching parameters are given in Table 2. These 
parameters were adjusted until numerical artifacts were eliminated from flow monitors such 
as the length of the recirculation region behind the injector (discussed in Section 4.3) and 
the position of the recompression shock (discussed in Section 4.4). Note that, for a = 0, the 
nodes are clustered close to the z a end of the zone and, for a = 0.5, the nodes are clustered 
close to both boundaries. 

For the long injector cases, the y-coordinates were generated in a similar manner. How- 
ever, the x-coordinates were generated using a compression function due to Thomas et al. 
[18]. In each zone there were N + 1 nodes (including end points) specified by Equation 1 


and 


V = Vo 


sinh(0r) — A?) 
sinh{ [A,,) 


+ 1 


( 5 ) 


where 


A - = 2 ln 


1 + (e^ - l)vo 


( 6 ) 


1 + (e -0 - 1)770 

and 770 is the value of 77 at which maximum compression occurs. Details of the zonal bound 
axies and stretching parameters are given in Table 3. 


3.3. Initial Conditions and Inflow Conditions 

The initial conditions within the combustor were set to T = 3001^, P = 300Po and zero 
velocity so as to approximate the conditions in T4 experiments in which the test section of 
the facility is evacuated to approximately 2 Torr. 

The inflow conditions were set to approximate the flow in the T4 facility fitted with a 
Mach 5 contoured nozzle. The nozzle supply pressure was approximately P t = 52 MPa and 
total enthalpy was, H, = SAM J/ kg. Details of the experiment are recorded in [15]. 

The test flow conditions, both in the shock-reflection region and the nozzle, were com- 
puted as quasi-steady flows. The nozzle supply conditions (behind the reflected shock) were 
estimated with the FORTRAN program “ESTC” [19] which incorporated an equilibrium 
model for air. From the shock reflection conditions, the gas expanded adiabatically (and in 
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chemical equilibrium) to the measured nozzle supply pressure of 52 MPa. Using these condi- 
tions, the flow at the exit plane of the nozzle was estimated with the quasi-one-dimensional 
program “NENZF” [20] in which the test gas consisted of a mixture of the species iVj, N, 
O 2 , O , NO, NO + and e~ . The gas was assumed to be in chemical equilibrium at the nozzle 
throat but a finite rate chemistry model with 11 reactions was used in the expansion region 
of the nozzle. 

The “steady” conditions at the end of the nozzle (and inlet to the combustor) were spec- 
ified as free-stream velocity 3670ms -1 , a static pressure of 21.5 kPa, a static temperature of 
1165-fif, and oxygen and nitrogen mass fractions of 0.2314 and 0.7686 respectively. Although 
the NENZF program determined a species mixture that included finite amounts of 0 and 
NO, only a nonreacting mixture of O 2 and JV 2 was used in the SPARK-2D computations. 

The “transient” inflow was constructed from the nozzle supply pressure trace shown in 
Figure 3, assuming a quasi-steady flow and using the ESTC and NENZF codes to compute 
the nozzle exit conditions for a number of values of P,. The time evolution of P s was then 
approximated by the curve fit 

P. = 7peofc(0.215 + 4.540£ — 22.148* 3 ), 0 < i < 0.28, 

= FpeaJfeCexK- 0 - 775 ^ - 0.28))), 0.28 < t < 2.26, (7) 

where Ppe a k = 62 MPa, t is in ms and P t is in MPa. The conditions at the inlet to the 
scramjet combustor were then related to P 3 with the following equations 


Pin(kPa) 

= 0.435P, - 1.55, 

( 8 ) 

Ti n (K) 

= 525P in 0256 , 

( 9 ) 

U in (ms ~ l ) 

= 350T m °' 333 . 

(10) 


The histories of the inlet conditions axe shown in Figure 8. 
3.4. Numerical Damping 


Some difficulty was experienced in starting the solution with such low pressure and large 
incident-shock Mach number (which was approximately 10 for Case 1). The difficulty ap- 
peared to be caused by the interaction of the strong incident shock structure with the newly 
forming recompression shock. To get the calculation started, the CFL number was linearly 
varied from 0.1 to 0.8 over the first 18000 time steps. For Cases 2,3,4 and 5, the coefficient 
for the artificial viscosity terms was set to a constant value of -1.0 while, for Cases 1 and 6, it 
was set to -1.2. To faithfully capture heat transfer, temperature smearing was not included. 
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3.5. Normalization Procedure for Transient Flow 


The flow quantities, q(x , t), computed in the transient Cases 2, 3 & 5 were normalized with 
upstream reference values to produce equivalent quasi-steady values, q(x,t ). The particular 
transformation used, based on the “hypersonic equivalence principle” was 

q{x,t) 


q(x,t) = 


( 11 ) 


q{x 0, t — Tdelay)’ 

where the time delay, T de i ay , was given as the nominal transit time from the reference position 

X ~~~ «Eq 


T delay — 


Vi, 


( 12 ) 


Essentially this means that individual parcels of fluid were followed downstream through the 
combustor and measured quantities were referenced to the parcels’ initial states. (See, e.g., 
Section 4.8 of [4].) 


4. RESULTS AND DISCUSSION 

In previous studies of impulsively started flows [21] [22] [23] [24] one principal result 
has been the time taken for the flow to approach steady state. This establishment time 
T e is usually combined with a characteristic velocity, U c , and a length, L c , to form the 
dimensionless parameter 

G=^, (13) 


U c r e 

L c ’ 

which represents the ratio of the flow establishment time to the time for flow to proceed 
through the domain. For the flow features in the generic scramjet, such as boundary layers 
and recirculation regions, G will have different values. For the impulsively started flow over 
a flat plate, where the boundary layer is laminar and L c is the distance from the start of 
the plate, G is approximately 3 [21]. When the boundary layer is turbulent G ~ 2. For the 
wake of a sphere, with L c as the length of the recirculation region, Holden [22] gives G — 30, 
based on the pressure measurements and G ~ 70 based on heat transfer measurements. 

In this section, six numerical simulations of the transient flow in a generic scramjet com- 
bustor will be examined with emphasis on the approach to steady state (or quasi-steady state 
as determined by the hypersonic equivalence principle). The following parameters will be 
used to determine flow establishment times: the boundary layer thicknesses at the end of the 
injector strut, the length of the recirculation region, the location of the recompression shock 
and, the pressure, heat flux and shear stress on the y — 0.02357m wall. The computation 
of the flow development for each of the cases was continued to a time of approximately lmj 
which is typical of the flow time in the current generation of pulsed-flow facilities. 
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4.1. Overview of the Solutions 


Figures 9 and 10 show the pressure fields for the short injector geometry with steady 
inflow (Case 1) and transient inflow (Case 2) respectively. Contour spacing is 2 kPa, These 
figures highlight the similarity between the steady inflow and the transient inflow solutions. 
The pressure distributions along the y = 0 and y = 0.02357m walls are shown on the right 
side of each figure. Note that the stretching of the y-coordinate (by a factor of 5) changes 
the appearance of the shock and Mach angles. 

At the start of the duct, there are a pair of weak shocks generated by the initially rapid 
growth of the boundary layer on the duct walls. Although not evident in the pressure field 
plots, boundary layers develop along the upper wall and the injector surface. The boundary 
layer on the injector separates at the step (i = 0.114m) and reattaches to the y = 0 
boundary some distance downstream. The expansion from the corner of the injector, and 
the recompression shock formed near the end of the recirculation zone are clearly visible for 
both cases. Subsequent reflections of the recompression shock from the duct walls can be 
seen for t > 0.2ms. 

The initial shock structure sweeping through the duct for t < 0.16ms is composed of an 
incident shock propagating downstream through the quiescent gas, a contact discontinuity 
separating the initial duct gas and the incoming test gas, and an upstream-facing shock 
decelerating the incoming test gas. Both the second-order and fourth-order MacCormack 
schemes experienced some difficulty with this strong shock structure and numerical oscilla- 
tions are evident in the solutions for t < 0.20ms. However, it is expected that the solutions 
are reliable at later times. 

The incident-shock position was identified as the point at which the wall pressure rose 
from (the initial value of) 300Pa to 800Pa. The position was recorded every 100/xs and the 
data points were then fitted with the straight line 

X = v(t - to). (14) 

The values of v and t 0 ar e given in Table 4. In all cases, the shock speed was virtually 
constant throughout the flow domain. However, the transient inflow cases (2,3 and 5) had 
a lower incident-shock speed than the steady inflow cases. For Cases 1 and 6, the position 
of the upstream-facing shock was identified as the point where the wall pressure crossed a 
35 kPa threshold. The velocities for this shock are also given in Table 4. 

As shown by the pressure contour plots and the wall pressure plots in Figures 9 and 
10, the overall flow features of both steady inflow and transient inflow cases are similar. 
Although the starting shock structure took longer to exit the duct for the transient inflow 
(Case 2), the major steady-state features such as the boundary- layer interaction shocks, the 
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expansion fan and the recompression shock are in virtually the same places. The principal 
difference seems to be the average level of the duct pressure which can be seen relaxing in 
the wall pressure plots of Figure 10. 

The effect of including the complete injector strut in the computational model can be 
seen in Figure 11, which shows a few frames from Cases 3 and 5. The frames are placed 
back to back and the times have been chosen to show approximately the same stages of 
development. However, the instantaneous inflow conditions are not exactly the same in each 
pair of frames. The major differences are the strong shock and expansion propagating from 
the front of the long injector. Also, the shock propagating from the leading edge of the 
y = 0.02357m wall is stronger than the boundary-layer interaction shocks seen in the short 
injector simulations. After reflecting from the injector strut and then the duct wall, this 
leading-edge shock interacts with the expansion fan propagating from the base of the strut. 

Since there is only a small difference between the long-injector and short-injector simula- 
tions, the following discussion will concentrate on the details of the short-injector simulations. 

4.2. Injector- Strut Boundary Layer 

For the short injector simulations, the boundary layer along the strut surface closely 
approximates the ideal flat-plate boundary layer, the main difference being the boundary- 
layer interaction shock striking the strut surface at x ~ 0.09m. This weak shock did not seem 
to affect the flow profiles within the boundary layer at the end of the strut (a; = 0.114m) 
where the following thicknesses were calculated: 


• total boundary layer thickness 5; the y - distance from the strut surface to 0.99 of the 
local free-stream velocity u e , 


• the displacement thickness 

*1 

1 

r—4 

II 

So 



(15) 

• the momentum thickness 




II 

£ £ 
1 

(16) 


The temporal development of 8, 8* and 0 is shown in Figure 12. For the steady inflow, 
the thicknesses approach to within 2% of steady state values by t = 0.23 — 0.26ms. If the 
approach to steady state begins with the passage of the upstream-facing shock, these times 
give 7 which is roughly double the value expected for a laminar boundary layer in shock- 
tube flow [21]. This result seems to indicate that the incident shock structure, including the 
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contact surface and the upstream facing shock, interferes with the establishment processes 
and delays the approach to steady state. 

The steady value of 5 — 1.49mm is close to the value of 1.45mm obtained with the 
“reference temperature” relations of Eckert [25] and Equation (7.53) in White [26]. This 
agreement provides some confidence that the viscous effects have been adequately resolved. 
Also, the flat plate transition data from He &; Morgan [27] indicate that transition to a 
turbulent boundary layer would be expected after a distance of approximately 0.25m. 

For the transient inflow cases, the approach to steady state is qualitatively different and, 
at late times, there is a slight growth in each of the thicknesses as a consequence of the 
relaxation of the free-stream pressure (and hence density). 

4.3. Recirculation Region 

The form of the velocity field in the recirculation region at the base of the strut appeared 
to be insensitive to the transient inflow. Figure 13 shows the velocity field in the recirculation 
region for Cases 2 and 3. The symmetry plane boundary condition in Case 3 allows a stronger 
vortex to form closer to the strut base. The pressures along the y = 0 line in the recirculation 
region are shown in Figure 14 at selected times. 

The length of the recirculation region, as measured by the position of the stagnation point 
near the y = 0 boundary, is shown in Figure 15. At late times, both the steady inflow and the 
relaxing inflow traces reach steady values. The steady state lengths are L r = 9.5, 8.3, 10.1 
and 10.8mm for Cases 1,2,3 and 6 respectively, which means that (a) the symmetry plane 
cases have longer recirculation zones than the solid wall cases and (b) the transient inflow 
cases have slightly shorter lengths than their steady inflow counterparts. Based on L r for 
steady inflow, the base flow establishes in a dimensionless time G — 146 ( t = 0.419ms) for 
the solid wall (Case 1) and G ~ 83 (t = 0.283ms) for the symmetry plane (Case 6). 

Also shown in Figure 15 is the length of the recirculation region along the line y = 
0.0025m which is half way up the base of the strut. This length was determined by computing 
the vertical mass flux 

TThy = ^ pUydX y (H) 

and finding the x-location where m v became zero. This measure showed the same trends 
as the stagnation point measurement except that G ~ 64 (< = 0.229ms) for the symmetry 
plane (Case 6). 

4.4. Recompression Shock 

Figure 16 shows the pressure distributions along the y = 0 plane for Cases 1 and 2. 
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The recompression shock is rather smeared and appears as the gradual pressure rise from 
x ~ 0.12m to i — 0.13m. Observe that the pressure distributions for the transient inflow 
(Case 2) collapse onto the one curve when they are normalized using the procedure discussed 
in Section 3.5. 

The shock position along a line through the midpoint of the strut base (y = 0.0025m) is 
plotted against time in Figure 17. The shock position was taken to be the grid location at 
which the pressure rises to 25% of the inlet pressure. Based on the length of the recirculation 
zone the nondimensional establishment time G ~ 140 ( t = 0.402ms) and G ~ 80 (t = 
0.273ms) for Cases 1 and 6 respectively. 


4.5. Wall Loads 


The shear stress and heat flux were monitored at selected points along the y = 0.02357 m 
wall. The shear stress coefficient is given by 


(i 

C f = ^ 


<9 tig 

a v wall 


and the heat flux by the Stanton number 


C h = 


%Pi n u in 2 ’ 


8T 

d JL mil 


Pin u in C p (T aw -T w y 


(18) 


(19) 


where the normal derivatives of the velocity and temperature were approximated by differ- 
ences at the first node off the wall. The viscosity was evaluated at the wall temperature using 
Sutherland’s formula and pi n and u in are inlet conditions computed with the time delay spec- 
ified in Equation (12). Thermal conductivity k was obtained assuming C p = 1004 J/ kg / K 
and Prandtl number Pr = pC p /k — 0.71. The adiabatic wall temperature was obtained 
from 


T — T 

-*■ aw — -* m 


l + VPr 


( 7 - 1 ) 


Mi, 


( 20 ) 


The development of these parameters at x — 0.114m and x = 0.5m is shown in Figures 18 
and 19 respectively. Results for the steady inflow cases exhibit an approach to steady state 
at both locations. The normalization of the transient inflow measurements is not entirely 
successful as the values of both Cf and Ch settle to a constant rate of growth after some 
start-up period. 

Normalized pressure histories at x — 0.114m and 0.50m wall are shown in Figures 18(c) 
and 19(c) for completeness. There is only a slight variation in the traces for the two transient 
inflow cases indicating that the normalization works well for pressures. The slight variation 
may be caused by the use of the constant value U c = 3670ms -1 in Equation (12). 
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4.6. (x,t)-Diagrams 


The time for wall pressure, shear stress and heat transfer to reach steady state at a number 
of x-locations along the y = 0.02357m wall is shown in Figures 20, 21 and 22 respectively. 
An average of each quantity between times t = 0.7 ms and 0.8ms was used as the steady 
value and the settling time t e determined such that the flow quantity remains within 2% of 
the steady value for t e < t < 0.7 ms. The establishment time is then given as r e = i e — t u f, 
where t u f t is the arrival time of the upstream- facing shock. Only the steady inflow cases (1 
and 6) are considered, since a consistent establishment criterion was not available for the 
transient inflow cases. 

The wall pressure generally settles more quickly than shear stress or heat transfer and 
has a dimensionless establishment time of G = T e U c jx ~ 0.9 for x > 0.2m. However, the 
establishment is significantly delayed where the relatively weak boundary-layer-interaction 
shocks strike the wall at x ~ 0.09m and 0.18m. The pressure measurement does not settle 
at all where the relatively strong recompression shock strikes the wall at x ~ 0.3m. 

The establishment times for the shear stress are nearly the same as those for the heat 
transfer and, for x > 0.4m, can be approximated by G ~ 3. A similar delay after the passage 
of the starting pulse was observed in the experimental measurements of East et al. [24]. For 
x < 0.3, the delay beyond the flat plate value of G ~ 3 is significant. 

5. CONCLUSION 

In order to study the establishment of the major flow features in a generic scramjet 
combustor, several numerical simulations were analysed. Of particular interest was the 
effect of a transient relaxing inflow condition as found in a reflected-shock tunnel which has 
been operated in an undertailored mode. 

From the numerical point of view, the simulations were particularly demanding and much 
of the interesting phenomena (such as turbulent mixing and chemistry) was omitted from 
these simulations. These issues will be addressed in a future study. The finite-difference 
scheme generally performed well but encountered difficulty with the very strong shock inter- 
actions and tended to produce oscillatory solutions at early times. With an incident-shock 
Mach number of approximately 10, these oscillations were not unexpected. An “upwinding” 
scheme (available in another version of the code) may better cope with the strong shocks. 

From the experimental point of view, it appears that the correlations available in the 
literature are adequate for predicting the establishment times for the flow features of the 
model if the inflow conditions are steady and allowance is made for the passage of the starting 
pulse [24]. For the relaxing flow, there are perturbations outside the 2% error band for a 
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significantly longer time but these perturbations are probably smaller than the measurement 
uncertainty in the shock tunnel experiment. Hence, relaxing test flow with a duration of the 
order of a few hundred microseconds is probably adequate for establishing flow in moderate 
sized scramjet models if the recirculating zones are not too large. The establishment time 
for turbulent jet mixing in the presence of combustion needs to be examined. 

The effect of the relaxing inflow conditions of the undertailored reflected-shock tunnel 
seem not to affect the overall pressure fields within the model nor the length of the recircu- 
lation region to any great extent. However, the relaxing flow does affect features such as the 
heat transfer and wall shear stress. It may also be expected to affect the chemical reaction 
rates which will scale as pL c and p 2 L c for two-body and three-body reactions respectively. 
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Table 1: Cases Considered in this Study 


Case 

Geometry 

Inflow for / > 0 

Spatial Differencing 

injector 

y = 0 boundary 

steady/transient 

1 

short 

cold wall 

steady 

2nd order 

2 

short 

cold wall 

transient 

2nd order 

3 

short 

symmetry plane 

transient 

2nd order 

4 

long 

symmetry plane 

steady 

4th order 

5 

long 

symmetry plane 

transient 

4th order 

6 

short 

symmetry plane 

steady 

2nd order 


Table 2: Grid Generation Parameters for Cases 1, 2, 3, and 6 


(a) y-grid 



Va 

(m) 

Vb 

(m) 

a 

P 

N 

Aj/min 

(m) 

■ 

0.005 

m 

0.5 


34 

50 x 10- 6 

H 

0.02357 

IB 

0.5 

DEES3 

51 

50 x 10-® 


(b) z-grid 



(m) 

Zb 

(m) 

■ 


N 

AXjnin 

(m) 

■ 

0.1 

0.0 

n 

1.153 


0.001 

■x 

0.114 

0.1 

M 

1.1495 

22 

0.00025 

n 

0.114 

0.124 

equal spacing 


0.00025 

n 

0.124 

0.168 

0.0 

1.042 

30 

0.00025 

B 

0.50 

0.168 ! 

equal spacing 

111 

0.00299 
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Table 3: Grid Generation Parameters for Cases 4 and 5 


(a) j/-grid for x > —0.10964m 


Zone 

3/a 

(m) 

yb 

(m) 

a 

p 

N 

^2/min 

(m) 

1 

0.02357 


1 

1 

H 

51 x 10- 6 

2 

0.0050 


i 

EH 

m 

53 x 10“ 6 


(b) y-grid for -0.138m < x < -0.10964m and x = (x + 0.138)/0.02836 


Zone 

3/a 

(m) 


a 

p 

N 

Aj/m in 

(m) 

1 

0.02357 

0.0050 x 

i 

B 

70 

65 x lO" 6 

2 

0.0050 x 

-0.00639 (1 - x) 

i 

EH 

30 

67 x 10" 6 


(c) x-grid 


Zone 

®a 

(m) 

Xb 

(m) 

Vo 

p 

N 

^®min 

(m) 

1 

-0.0813 

-0.138 


1.5 

30 

0.00173 

2 

0.0813 

-0.0813 

'/r-M 

3.5 

110 

0.00093 

3 

0.1467 

0.0813 


7.0 

80 

0.00017 

4 

0.2044 

0.1467 

AXt 

Ax;_i 


19 

0.00277 

5 

0.50 

0.2044 

equal spacing 

89 

0.00332 
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Table 4: Incident Shock (is) Speed and Upstream-Facing Shock (ufs) Speed as defined by 

Equation (14). 


Case 

Vi, 

(m/ms) 



tufa 

(ms) 

1 

3.68 

0 

2.83 

0 

2 

2.96 

0 

- 

- 

3 

2.97 

0 

- 

- 

4 

3.72 

0.0350 

- 

- 

5 

2.97 

0.0466 

- 

- 

6 

3.71 

0 

2.83 

0 
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Figure 1: Free-piston driven reflected-shock tunnel. 

(a) Facility schematic. 

(b) Distance-time wave diagram after primary diaphragm rupture. 


(a) Pitot pressure, Ppjtot 

2.0 

1.6 

Ppitot- 1-2 
MPa -8 
.4 
0 

2.0 2.5 3.0 3.5 4.0 4.5 

t, ms 



t, ms 

Figure 2: Typical pressure histories for the T4 shock tunnel. 

(a) Pitot pressure at the exit plane of a Mach 5 nozzle. 

(b) Nozzle supply pressure. 
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(a) 


Driver tube 


Shock tube 


Acceleration tube 


Primary Secondary 
diaphragm diaphragm 


(b) 



Unsteady 
expansion 
Test 


region 

Contact 

Shock surface 


Figure 3: Expansion tube facility (a) and wave diagram (b). 



Figure 4: Typical static pressure in the HYPULSE facility. 
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Test flow 



Figure 5: A generic scramjet combustor with centrally located fuel injector strut. 


(a) Short injector configuration 



Cold wall 


(b) Long injector configuration 

C H 

H — 0.138m — *+•— 0.114m— H 

h 0.5m 


Figure 6: Combustor geometry used in the computations. 
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Figure 7: Grids used for the calculations. Note that the y-scale is magnified by a factor of 
five. 




Static pressure 




time (ms) 


Figure 8: Transient inlet conditions as specified by equations (7) - (10). 













































t = .16ms 


t - .206ms 
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Case 3, transient with symmetry plane 


t = 0.1 ms 


0116 x, m °- r 
Case 2, transient with no-slip wall 


^ t = 0.1 ms 



0.4, 0.6, 0.8 


0.116 0.118 
x, m 


Figure 14: Pressure distributions along y = 0 near the base of the injector for Cases 2 and 
3. 


l p mm 


10 - 
5 - 



L , mm I 
2 



Figure 15: Length of the recirculation region. 

(a) L t identified by the stagnation point location. 

(b) length identified by a zero-crossing of the vertical mass flux along the line y = 0.0025m 











Figure 16: Pressures distribution at y = 0 downstream of the injector for Cases 1 and 2. 
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0.8 
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0.2 


0.4 


0.6 


Normalized pressure 


0.8 1.0 


t, ms 


Figure 18: Comparison of loads on the y = 0.02357m wall at x = 0.114m for steady and 
transient inflow. 
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Figure 19: Comparison of loads on the y = 0.02357m wall at x = 0.5m for steady and 
transient inflow. 








Figure 20: Establishment time for normalized pressure along the y = 0.02357m wall. 


IjO 
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0.3 0.4 0.5 

x, m 


coefficient along the y = 0.02357m wall. 
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